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Zusammenfassung 

We present a method for time series analysis of both, scalar and nonscalar 
time-delay systems. If the dynamics of the system investigated is governed 
by a time-delay induced instability, the method allows to determine the delay 
time. In a second step, the time-delay differential equation can be recovered 
from the time series. The method is a generalization of our recently proposed 
method suitable for time series analysis of scalar time-delay systems. The 
dynamics is not required to be settled on its attractor, which also makes 
transient motion accessible to the analysis. If the motion actually takes place 
on a chaotic attractor, the applicability of the method does not depend on 
the dimensionality of the chaotic attractor - one main advantage over all 
time series analysis methods known until now. For demonstration, we analyze 
time series, which are obtained with the help of the numerical integration of 
a two-dimensional time-delay differential equation. After having determined 
the delay time, we recover the nonscalar time-delay differential equation from 
the time series, in agreement with the 'original' time-delay equation. Finally, 
possible applications of our analysis method in such different fields as medicine, 
hydrodynamics, laser physics, and chemistry are discussed. 

P.A.C.S.; 05.45. -fb 



1 Introduction 

Time-delay differential equations have been widely proposed to account 
for the observed oscillatory, chaotic or hyperchaotic motion of dynamical 
systems. Most of the attention has been devoted to scalar time-delayed 
models |T]] - 0- Since the pioneering work of Farmer it is well 
established that scalar time-delay differential equations are able to 
exhibit high- dimensional chaotic attractors with many positive Lyapunov 
exponents and, therefore, are prominent examples to illustrate the chaotic 
hierarchy JTOf- Since then, scalar time-delay equations, especially the 
well-studied Mackey-Glass system IQ , have been used as model systems 
to produce high- dimensional chaotic time series. In the case of small 
delay times, the resulting low-dimensional chaotic dynamics is accessible 
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to time series analysis with the help of well-established methods |12 |- 
pif . For example, the fractal dimension and the Lyapunov exponents 
of the chaotic attractors can be estimated. In the case of large delay 
times, where the dynamics is high-dimensional chaotic, these methods 
run into severe problems [ P^ . 

A first step towards the time series analysis of time-delay systems 
has been done by Fowler and Kember who showed how 'smart 
embeddings' can indicate the presence of an underlying scalar time- 
delay system for any delay time. Later, we introduced a time series 
analysis method in order to verify the existence of an underlying time- 
delay system [O]-|20|. If the dynamics is governed by a scalar time-delay 
differential equation, we have shown that the delay time and the time- 
delay differential equation can uniquely be recovered from the time 
series. The method does not put any restriction on the dimensionality 
of the dynamics analyzed, opening up a door towards the time series 
analysis of high-dimensional chaotic motion in time-delay systems. Furthermore, 
the method does not require the motion to be settled on its attractor. 
The method turned out to be practically insensitive against additional 
noise, hence providing a well-suited tool for the analysis of experimental 
time series. So far, we have successfully applied the method to time 
series taken from an electronic oscillator [|T^ and a computer experiment 



|T9| . Although some of the time-delay models investigated are indeed 
scalar ones, time-delay systems are, in general, nonscalar. Nonscalar 
time-delay differential equations have been proposed in such different 
fields as arrays of coupled oscillators [^, laser physics PB[]-P^, 
physiology -[^; hydrodynamics |^|, and chemistry to account 
for the observed unstable and chaotic dynamical behavior. Additionally, 
several models with multiple delay times have been investigated pl] - 



33| . In this paper, we present a generalization of the time series analysis 



method proposed in ||T^-||T9[. Herewith, we are able to identify nonscalar 
time-delay systems and, therefore, verify the existence of an underlying 
time-delay induced instability. The method allows the recovery of the 
nonscalar time-delay differential equation from the time series. This 
article is organized as follows. In Section 2, the basic idea of the time 
series analysis method is presented. Also, adequate measures to determine 
the delay time from the time series are discussed. In Section 3, we 
illustrate the method by applying it to a two-dimensional time-delay 
system, the trajectories of which are computed numerically. We show 
that, while the delay time can still be determined even with a scalar 
ansatz, the existence of an underlying time-delay system cannot be 
verified. The latter is accomplished with the help of a nonscalar ansatz, 
which additionally allows to recover the time-delay equation from the 
time series. Finally, possible applications to experimental systems are 
discussed in Section 4. 
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2 The Basic Idea of the Analysis 



We consider an A^-dimensional time-delay differential equation 

Ut) = hiyo{t),yroit)), (1) 

with the initial condition 

yo{t) = yi{t), -To < t < 0. (2) 

The time evolution of yo at time t does not only depend on its present 
state ^0(^)5 it also depends on a state in the past, yr^ (t), which introduces 
nonlocal correlations in time. The state of the system is uniquely defined 
by functions on an interval of length tq. Therefore, the phase space 
of system (|I|) is C^, where C^-o is the space of continuous functions on 
the interval [— tq, 0] and the phase space has to be considered as being 
infinite dimensional. 

The trajectory in the infinite dimensional phase space 3^(t, t) G 
can be recovered from its projection yo{t) without loss of information 

y{T,t) = yo{t-T), -To<r<0. (3) 

We emphasize that the construction (^) of the trajectory in phase space 
is exact and can be accomplished for all values of the control parameters 
of the time-delay equation. The N time series ^o(^) encompass the 
complete information about the dynamics of the system in the infinite 
dimensional phase space C^, no matter which dynamical state is realized 
by the time-delay system. Equation (Q) holds for transient motions as 
well as for motions on chaotic attractors of arbitrary dimension. In 
general, it is expected that the dynamics of an infinite dimensional 
system cannot be reconstructed from a finite number of time series. 
For instance, to construct the trajectory in the phase space of spatial 
systems, the dynamics of which is governed by nonlinear partial differential 
equations, an increasing number of time series with increasing complexity 
of the dynamics is required. 

Additionally, the time derivative y{T,t) G of the trajectory in 
phase space can be estimated in principle 

3>(r,t) = y'o(t-r), -ro<r<0, (4) 

though, in practice, the estimation of time derivatives can be the source 
of severe errors. 

In the case of ordinary as well as time-delay differential equations, 
the time derivative of the trajectory in phase space is functionally 
related to the trajectory in phase space via the time-evolution equation. 
It is the specific property of time-delay systems, though, that only a 
restricted number of coordinates are correlated via the time-evolution 
equation (|lD, namely, the 2N coordinates {yo{t) , yro{t)) taken from the 
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trajectory y{r,t) in phase space with the N coordinates yo(t) taken 
from the time derivative 3^(t, t) of the trajectory. 

The basic idea of the time series analysis method presented is to test 
whether a time-delay equation (|l|) can be constructed, the solution of 
which is given by the observed time series. In order to accomplish this 
task, in general, one has to construct a phase space with the help of the 
observed time series. In infinite dimensional systems, one would expect 
that an increasing number of time series with increasing complexity of 
the dynamics is required. We have argued in the preceding section that 
this is not true for time-delay systems, but the observation of time 
series is sufficient to construct the trajectory in the infinite dimensional 
phase space of an A^-dimensional time-delay system, no matter which 
dynamical state is realized. Secondly, it is the specific property of 
time-delay systems that only a restricted number of coordinates are 
correlated via the time-evolution equation namely, the 2N coordinates 
{yo{t) , yrg{t)) taken from the phase space and one of its time derivatives 
Voit). Therefore, it is not necessary to analyze the dynamics in the 
infinite dimensional phase space, to verify the existence of an underlying 
time-delay system. It is sufficient to show the existence of a functional 
relationship (p. To this end, we analyze the dynamics in a 3A^-dimensional 
space, which is spanned by the coordinates {yo, yr^, yo). The dynamics of 
a time-delay system in the 3A^-dimensional space is restricted to a 2N- 
dimensional hypersurface, which is given by the time-evolution equation 
([T|). The hypersurface (|I|) is defined that to any value (^0(^)5 yroif)) there 
is a unique value of the time derivative yo{t) for all times t. 

The idea of the time series analysis method presented in this article 
is to test the existence of such a hypersurface for a given time series. 
Starting with N scalar time series, which have been taken from the 
system to be investigated, we hypothesize, at first, that the dynamics 
is governed by an A^-dimensional time-delay system 

Voit) = hr{yo{t),yT{t)), (5) 

with an unknown function hr and an unknown delay time r, both of 
which will be determined in the subsequent analysis, if the ansatz (^) 
turns out to be successful. Then, we take the values of {yo,yr) and yo 
from the time series and analyze its dynamics in the 3A^-dimensional 
space, which is spanned by the coordinates {yo, yr, yo). If the coordinates 
of the trajectory {yo{t) , yrit) , 2/0 (^)) functionally correlated via equation 
(|^), the hypothesis that the system is governed by a time-delay equation 
with the delay time r has been verified. If the projected trajectory 
{yo, yr, yo) does not fulfill condition (^, the hypothesis has to be rejected. 
This is a unique criterion to determine the delay time from the time 
series. Additionally, the time-delay differential equation can be constructed 
by analyzing the functional relationship (^), which exactly gives the 
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function h^. Obviously, the only requirement remains that the system 
dynamics fulfills the time-evolution equation which even is true for 
all kinds of transient motion as well as for the motion on chaotic or 
hyperchaotic attractors of arbitrary dimension. Therefore, the method 
permits to analyze high-dimensional chaotic dynamics of time-retarded 
systems, which is not accessible to the fractal dimension analysis. If 
the system possesses several coexisting attractors 0, the method is 
applicable to the dynamics on every attr actor. 

We have argued above that the existence of an underlying time-delay 
system can be verified by proving the existence of a 2A^-dimensional 
hypersurface hr in the form in the 3 A^- dimensional space, which is 
spanned by the coordinates {yo,yr,yo)- Therefore, it is crucial to apply 
adequate measures, which enable us to identify such a hypersurface by 
analyzing the time series. All measures proposed so far |]l6l - [|19| solely 
rely on the fact that, if the trajectory {yo{t),yr(t),yo(t)) is correlated 
via equation (|]), the dimensionality of the trajectory is reduced. To 
our knowledge, this has been realized the first time by Fowler and 
Kember [|T6[, who analyzed the dynamics of the Mackey-Glass equation. 
They applied an embedding of the time series in a three-dimensional 
space with two time-delayed coordinates. The delay time of the first 
coordinate has been chosen to be small. The delay time of the second 
coordinate was taken as variable. The authors of [^] stated that if 
the delay time of the second coordinate equals the delay time of the 
time-delay system, the trajectory lies 'close to a surface'. Fowler and 
Kember applied a singular value fraction to detect the decrease in 
dimensionality. As has been correctly mentioned by them, the singular 
value fraction is not a good tool, if the surface is folded. The latter 
must be considered as the general case. 

Recently [0]"I13' have proposed a time series analysis method 
for scalar time-delay systems, only. There, we showed that the trajectory 
in the (t/q, ?/to! J'o)-space is restricted to a two-dimensional surface. The 
reduction in dimension was detected by intersecting the projected trajectory 
with a surface k{yo,yr,yo) = 0. The intersection points must be on a 
curve, if the projected trajectory is correlated via equation (H). We 
have detected such a behavior by ordering the intersection points with 
respect to one coordinate and drawing a polygon line, which connects 
all intersection points. The length of the polygon line has been taken 
as a measure for the alignment of the intersection points. This simple 
method correctly determined the delay time of scalar time-delay systems, 
but there is no straightforward generalization for nonscalar time-delay 
systems. 

For this reason, we apply another method, the basic idea of which is 
the following: If the trajectory of an A^-dimensional time-delay system 
in the {yo, yVo, yo)"Space, where tq is assumed to be the correct value of 
the delay time, is correlated via equation (|^), the trajectory is restricted 
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to a hypersurface. Therefore, most parts of the {yo,yrg,yQ)-spa.ce are 
not visited by the trajectory. If the trajectory is viewed in any other 
space, for instance, if the value of the delay time is not chosen properly, 
the projected trajectory is expected to visit 'more' parts of the space. 
Therefore, we compute the filling factor of the projected trajectory by 
covering the {yo,yr,yo)-spa.ce with equally sized hypercubes. The 
filling factor is the number of hypercubes, which are visited by the 
projected trajectory, normalized to the total number of hypercubes, 
P^^ . The filling factor is computed under variation of r. The existence 
of an underlying time-delay induced instability induces a local minimum 
in the filling factor. 

Fowler and Kember |T^ already suggested a fractal dimension analysis 
of the projected trajectories. The fractal dimension analysis has some 
severe drawbacks: At first, a fractal analysis requires a large number of 
data points, because for the determination of a fractal dimension it is 
crucial to resolve the geometrical object under investigation on different 
'length scales'. Secondly, the fractal analysis is computationally intensive 
and, in practice, sensitive to additional noise. All those measures rely 
on the fact that the trajectory is restricted to a hypersurface, if the 
projected trajectory is correlated via equation (|^). It has been argued 
above that the dimension reduction is not a sufficient criterion for the 
verification of an underlying time-delay system. The existence of the 
functional relationship has to be shown separately. 



3 Time Series Analysis of a Two-Dimensional Time- 
Delay System - A Numerical Example 

The applicability of the method for time series analysis is demonstrated 
with the help of a computer experiment. A nonscalar time-delay system 
is integrated numerically. Details of the numerical integration are reported 
in the appendix. We will show that, although a scalar ansatz also leads 
to a local minimum in the filling factor, the scalar ansatz must be 
rejected, because it is not possible to find a surface, which is given by 
an equation of the form with the help of a scalar ansatz. In a second 
step, we will identify the system as a nonscalar time-delay system by 
verifying the existence of such a surface with the help of a nonscalar 
ansatz. Finally, the time-delay differential equation will be recovered 
from the time series. 

We consider the two-dimensional time-delay differential equation, 
which has been chosen to serve its demonstrational purpose best: 

u = -V + f{Uro), (6) 

V = g{u,v), 
with the initial condition: 

uit) = -ro<t<0. 
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v{t = 0) = Vi. 



The functions / and g are given by: 



g{u,v) 




(7) 
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Equation (|^) has some similarity with the Mackey- Glass system. The 
dependence of ii on the time-delayed value Uro is the same as it is the 
case in the Mackey-Glass system. But while the dependence of -u on m 
induces exponential relaxations in the Mackey-Glass system, in it 
is similar to a damped oscillator. There are two limits of the control 
parameter space (a, T, tq), where the dynamics of @ is well-known. For 
a ^ 0, system (P) reduces to a damped harmonic oscillator. For T — > 0, 
the time scale of v is much faster compared to the time scale of u. The 
variable v, then, adiabatically follows variable u and the dynamics of (P) 
resembles that of the Mackey-Glass system. Here, a and tq are chosen 
such that in the scalar limit (T —> 0) the dynamics is high-dimensional 
chaotic (a = 3,to = 20). T is varied from 0.10 to 1.90. Note that the 
system |^ can be transformed to a scalar integro-differential equation 
for the variable u by integrating the second equation with the help of 
the method of varying coefficients. 

We present three time series of u and v for different T in Fig. 1. In 
Fig. 1(a), it is clearly seen how the variable v, for T = 0.10, follows 
the variable u and the dynamics of u resembles that of a Mackey-Glass 
system. The values of {u, v) are positive for all times. We mention at this 
point that system (|]) is invariant under the transformation {u, v) — > 
{—u, —v). Therefore, there exists another attractor with negative values 
of {u,v). In Fig. 1(b) and Fig. 1(c), we observe that the two coexisting 
attractors are merged. Variable v no longer follows variable u, but it 
develops an independent dynamics. In these cases, the system reveals 
its nonscalar nature. 

Now, we analyze these time series with the help of a filling factor 
analysis. At first, we choose the scalar ansatz 



for the analysis of the time series u{t) with an unknown delay time 
T and an unknown function /i^.We analyze the time series in a three- 
dimensional space, which is spanned by the coordinates [u, Ur, ii) with 
a variable value of r. Then, the filling factor of the time series is 
determined under variation of r. 

The results of the filling factor analysis are presented in Fig. 2 for 
different values of T. The filling factor is minimal for small values of 
r as a result of local correlations in time. The filling factor increases 
for increasing r, which is a fingerprint of the chaotic nature of the 
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motion, and eventually reaches a maximal value. For r = Tq = 20.00, a 
local minimum of the filling factor is observed for all values of T. This 
decrease in the filling factor is due to the nonlocal correlations in time 
induced by the time delay. An additional local minimum appears at 
r = 2ro. For high enough values of T, though, other regularly spaced 
local minima in the filling factor appear, the period of which is equal 
to the oscillations of the underlying damped oscillator. In Fig. 2(b), a 
blow-up of the r-dependent filling factor in the vicinity of the delay 
time tq = 20.00 is shown. Clearly, the local minimum appears for all 
values of T considered here, but it is less pronounced for increasing 
T, because the character of the time-delay system is more and more 
becoming nonscalar, then. As has been emphasized above, the reduction 
in dimension is only a necessary, but not a sufficient condition for the 
verification of a scalar time-delay system. The existence of a surface in 
the form has to be checked as well. 

To this end, we apply an intersection of the time series u{t) with 
the plane u{f) = 0. The values of the coordinates = u{f) and 
u^.^ = Uro{f) are recorded. They have to be correlated according to the 
scalar ansatz (^: 



Therefore, the points have to lie on a curve, if the dynamics 

is governed by a scalar time-delay equation. The intersection points 
(m*, ul-J for r = Tq = 20.00 and different values of T are shown in Fig. 
3. For small values of T, the time series v{t) follows the time series 
u{t). In this case, the dynamics of variable u is close to the dynamics 
of the Mackey-Glass system. Such a behavior can be seen in Fig. 3(a), 
where the intersection yields a geometrical object, which is close to 
being a one-dimensional curve. The inset of Fig. 3(a) shows a blow- 
up of the intersection points. Clearly, the alignment of the intersection 
points is not perfect, as a result of the nonscalar nature of system (H). 
Nevertheless, the scalar ansatz (^ would be a good approximation for 
small values of T. For higher values of T, system (^ reveals its nonscalar 
nature. Therefore, in Fig. 3(b) and Fig. 3(c), the distribution of the 
intersection points becomes cloudy. That means, a smooth functional 
relationship (^ cannot be found in these cases and the scalar ansatz 
has to be rejected. Nevertheless, in the considered example of a two- 
dimensional time-delay system the filling factor analysis was successful 
in the framework of a scalar ansatz, in the sense that the r-dependent 
filling factor showed local minima for the correct values of the delay 
time. 

Now, we will analyze the time series for higher values of T with the 
help of a nonscalar ansatz. The general two-dimensional ansatz is 



with an unknown delay time r and two unknown functions hr^i,hr,2- 




(10) 





(12) 
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The analysis has to be conducted in a six-dimensional space, in which 
the dynamics is restricted to a four- dimensional hypersurface. The 
hypersurface is given by the functions /i,. i, /ir,2, which can be determined 
with the help of adequate fitting procedures, for instance, a least- 
squares-fit in the framework of a presupposed model. In this article, 
we chose a more restrictive ansatz for demonstrational purposes, 

U = -V + fr{Ur), (13) 

V = -gr{u,v). (14) 

The delay time r and the functions fr and gr are yet unknown and 
will be determined in the following. At first, we perform a filling factor 
analysis in the same spirit as has been done in the scalar case. But now 
the two time series {u{t),v{t)) have to be projected to a six-dimensional 
space which is spanned by the coordinates {u,Ur,u,v,Vr,v)- The six- 
dimensional space is covered with equally sized hypercubes and the 
number of hypercubes which have been visited by the trajectory is 
counted under variation of r. 

The results are presented in Fig. 4 for different values of T, namely 
T = 0.10,0.60,1.90. The minimum in the filling factor for t = Tq = 
20.00 is well detected for all values of T. In Fig. 4(b), we show a blow- 
up of the r-dependent filling factor in the vicinity of the delay time Tq 
of system (|^). Clearly, the local minimum for r = tq is detected for all 
values of T. 

As has been argued above, the existence of the functional relationship 
d^) has to be shown, in order to verify the underlying time-delay induced 
instability. The special form of ansatz (O)-(n) together with the time- 



evolution equations (^ allows for a convenient way of proving the 
existence of the function (^). We emphasize, though, that, in general, 
it is expected to be more troublesome. We apply an intersection with 
the help of the condition v{f) =0. If the nonscalar ansatz (|T3p-(0) is 
successful, the values -it* = iiit'') and u].^ = Ur^if) have to be correlated 
via 

= /.«). (15) 

Plotting versus u^.^ as is shown in Fig. 5(a), the existence of the 
smooth function fr is verified. We compare the reconstructed function 
fr (open circles) with the function / (line) of the time-delay differential 
equation @. The coincidence is good and the {u\ u^J-plot can be used 
to recover the function / from the time series. We emphasize that 
no parameter has been adjusted to compare the function / with its 
recovery fr in Fig. 5(a). 

In the next step, the functional relationship between and is 
investigated. We use the same intersection condition v{f) = as above. 
According to the nonscalar ansatz (|13|)-(|14|), the coordinates of the 
intersection points are correlated via 

v' = -griu\0). (16) 
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Plotting versus yields the function gr{u\0). We compare the 
functions g{u\0) and gr{u\0) in Fig. 5(b). Again, the correspondence 
is good and the {v\ ti*)-plot can be used to recover the function gr{u\ 0) 
from the time series. The recovery of the function gr{0, v^) is to be done 
in the same spirit and is not shown here. Obviously, the existence of the 
functions fr and gr has been proven and a two-dimensional time-delay 
system has been identified by analyzing the time series. 



4 Applications 

Finally, we would like to discuss possible applications of the present 
method for time series analysis. We emphasize that our method does not 
have the restrictions which are inevitable to the embedding techniques 
necessary for the determination of fractal dimensions of chaotic attractors 
in phase space. The analysis is not restricted to a low-dimensional 
chaotic motion. Transients can be analyzed as well. The method is 
not sensitive to additional noise. Furthermore, it has been shown that 
the analysis can be performed with a comparably small number of data 
points [^]. Apparently, we have a well-suited tool for the analysis of 
experimental time series. If the dynamics of the system to be investigated 
is governed by a time-delay induced instability, the method allows for 
the determination of the delay time and a recovery of the time-delay 
differential equation. Therefore, it is possible to compare time series of 
experimental systems with proposed model equations in detail. System 
parameters can be extracted from the time series analysis, which might 
be not accessible otherwise. 

We speculate that the analysis can be particularly useful in such 
fields as medicine and biology, where noninvasive techniques are of 
great importance for obvious reasons. In several experiments on human 
subjects, which were exposed to time delays of some sort, a qualitative 
change of the observed dynamics has been verified [||, ||, |2^, which 
could be correlated to well-established pathologies. If the observed 
dynamics is, indeed, induced by a time delay, as proposed, we expect 
our analysis method to be successful. On the other hand, we propose 
to check the validity of the Mackey-Glass system by analyzing suitable 
time series. The analysis might improve the understanding of the experiments 
and possibly allow to determine important system parameters and serve 
as a new diagnostic tool. 



In a recent article, Villermaux deals with the low-frequency 
oscillations of the velocity field in the 'hard turbulence' regime in a 
closed convection box. The author proposes a two-dimensional time- 
delay system to describe the dynamics of the disturbances. Our method 
has the potential to test the validity of the model by analyzing the 
experimental time series, which possibly can lead to a better understanding 
of boundary instabilities. 

Another important class of a prototype model for chaotic dynamics 



10 



are laser systems. No wonder that time-delay models have also been 
investigated in laser physics ([^ ^ and references therein). 

It is the advantage of laser systems with a time delay that they allow 
precise measurements of time series p5| , ^ . We find them particularly 
suitable for our analysis, because we expect the time delay of laser 
systems to be practically discrete, compared to other experimental 
systems. The first steps towards the identification of high- dimensional 



chaotic dynamics of a time-delayed laser system has been taken [20 



5 Concluding Remarks 

In conclusion, we have presented a generalization of a recently proposed 
method for recovering the time-evolution equation of scalar time-delay 
systems by analyzing the time series. The method is generalized in the 
way that it can be applied to nonscalar time-delay systems. 

We have shown that an A^-dimensional time-delay system can be 
identified with the help of time series. The analysis has not to 
be conducted in the infinite dimensional phase space, instead it is 
sufficient to analyze the dynamics in a 3iV-dimensional space, in which 
the dynamics has to be restricted to a 2A^-dimensional hypersurface. 
This finding gives us a unique criterion to determine the delay time of a 
nonscalar time-delay system by analyzing the time series. Additionally, 
the time-delay differential equation of dimension can be recovered. 

We emphasize that we only require the motion to obey the time- 
delay differential equation. The motion is not required to be located 
in certain parts of phase space. If the dynamical system possesses 
coexisting attractors, the method can be applied to motions on every 
coexisting attractor. Moreover, the method is also applicable to transient 
motions. If the motion is on a chaotic or hyperchaotic attractor, the 
applicability of the analysis method does neither depend on the dimensionality 
nor on the number of positive Lyapunov exponents of the chaotic or 
hyperchaotic attractor. Therefore, we find that the present method 
might open up a door towards the time series analysis of high-dimensional 
chaotic motion in time-delay systems. 

We have shown the applicability of the method by analyzing time 
series which have been obtained with the help of numerically integrating 
a two-dimensional time-delay differential equation. The system investigated 
mimics a scalar time-delay system in a certain parameter range, where 
a scalar ansatz yields a good approximation. Under variation of a single 
control parameter, the dynamics increasingly reveals the nonscalar nature 
of the time-evolution equation. In this scalar ansatz is not 

sufficient, but, nevertheless, we have successfully analyzed the time 
series with the help of a nonscalar ansatz. The delay time has been 
determined from the time series. Finally, we recovered the two-dimensional 
time-delay differential equation by analyzing the time series. Possible 
applications for the analysis of dynamical systems in such different fields 
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as medicine, hydrodynamics, laser physics, and chemistry are discussed. 
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7 Appendix: Numerical Methods 



The results presented in this article have been computed with the help 
of a Runge-Kutta algorithm of fourth order. The fundamental time step 
was taken to be 0.01. The length of the memory has been 2,000 time 
steps. In all simulations presented, we have chosen homogeneous intial 
values u{t) = 0.50, v{t) = 0.90, -Tq < t < 0. 

We checked the validity of the numerical integration of system (P) 
by comparing the results of Runge-Kutta algorithms of different order, 
under variation of the time step of the integration, in order to check 
the validity of the results. Additionally, the computed time series were 
compared to analytical solutions, which are available in certain parameter 
ranges. 

For the filling factor analysis, we used time series with 1,000,000 
data points. The time derivatives were estimated by applying a local 
parabohc approximation. To compute the r-dependent filling factor, 
every tenth point of the time series was taken. Therefore, we conducted 
the filling factor analysis by analyzing 100, 000 data points in the three- 
dimensional or six-dimensional space, respectively. 
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Figure captions 

Fig. 1: Time series of the system obtained with the help of numerical 
integration (a = 3.00, tq = 20.00) for different values of T: (a) 
T = 0.10; (b) T = 0.60; (c) T = 1.90. 

Fig. 2: (a) Filling factor of the time series u{t) of the system (j^) with 
a = 3.00, To = 20.00 for different values of T (lower curve: T = 
0.10; middle curve: 0.60; upper curve: 1.90). We used 100, 000 data 
points for the filling factor analysis, which were taken out of a time 
series of 1, 000, 000 data points. In (b) a blow-up of the filling factor 
in the vicinity of the delay time r = tq is shown [squares : T = 
0.10; circles : T = 0.60; stars : T = 1.90). 

Fig. 3: Intersection points of the time series u{t) for a = 3.00, tq = 20.00 
and different values of T: (a) T = 0.10, (b) T = 0.60, (c) T = 1.90. 
The analysis has been conducted with a time series of 1,000,000 
data points. 

Fig. 4: (a) Filling factor of the two time series {u{t),v{t)) of the system (^) 

with a = 3.00, tq = 20.00 for different values of r(T = 0.10, 0.60, 1.90; 
the values of T are indicated in the figure). We used 100,000 
data points for the filling factor analysis, which were taken out 
of a time series of 1, 000, 000 data points. In (b) a blow-up of the 
filling factor in the vicinity of the delay time r = tq is shown 
[squares : T = 0.10; circles : T = 0.60; stars : T = 1.90). 

Fig. 5: Recovery of the time-delay differential equation from the time 
series: (a) Comparison of the data points [u\ u^^), which are shown 
as open circles, with the function / (line). For clarity, only 300 
data points [u\ u^.^) are shown, (b) Comparison of the data points 
(v*, M*), which are shown as open circles, with the function g[u^, 0) 
(line). For clarity, only 40 data points [v\u^) are shown. 
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